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We investigate the abundance of large-scale hot and cold spots in the WMAP-5 temperature 
maps and find considerable discrepancies compared to Gaussian simulations based on the ACDM 
best-fit model. Too few spots are present in the reliably observed cosmic microwave background 
(CMB) region, i.e., outside the foreground-contaminated parts excluded by the KQ75 mask. Even 
simulated maps created from the original WMAP-5 estimated multipoles contain more spots than 
visible in the measured CMB maps. A strong suppression of the lowest multipoles would lead 
to better agreement. The lack of spots is reflected in a low mean temperature fluctuation on 
scales of several degrees (4°— 8°), which is only shared by less than 1% (0.16%-0.62%) of Gaussian 
ACDM simulations. After removing the quadrupole, the probabilities change to 2.5%— 8.0%. This 
shows the importance of the anomalously low quadrupole for the statistical significance of the 
missing spots. We also analyze a possible violation of Gaussianity or statistical isotropy (spots 
are distributed differently outside and inside the masked region). 
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I. INTRODUCTION 

The precise measurement of anisotropics in the cosmic 
microwave background (CMB) has played a key role 
in amplifying our knowledge about the structure and 
evolution of the Universe. The best data available today 
is provided by the Wilkinson Microwave Anisotropy Probe 
(WMAP) satellite mission from five years of observation. 
Its results arc powerful enough to put various cosmological 
models to stringent tests. They helped establishing 
today's standard model of a spatially flat universe with 
Gaussian initial perturbations, possibly generated during 
an early inflationary epoch. According to the standard 
ACDM model, the present Universe is essentially made up 
from dark energy in the form of a cosmological constant 
A and cold dark matter (CDM). Under the assumptions 
of Gaussianity and statistical isotropy, all the information 
about the temperature fluctuations in the CMB are en- 
coded in the angular power spectrum Ci from a harmonic 
decomposition of the temperature field. A crucial result 
of the WMAP analysis therefore is an estimate of the 
multipol es Ce which i s in g ood agreement with the ACDM 
best fit ( Nolta et all I2009T ) except for the well-known dis- 



crep ancies of the low mult ipoles, especially the quadrupole 
C 2 (|Hinshaw et all . 120071) . Nonetheless, many issues are 



still under intense discussion. Repeatedl y, authors have 
claim ed to detect non-Gaus sia n sign als ( McEwen et all 



or statistical 



anisotropy (Bernui et al., 120061: Ic 


e Oliveira-Costa et al. 


2004; Eriksen et al., 2004 


Hansen et al, 2009; 


Hoftuft et al., 2009t Land and Maeuciio. 2005). Since 



the power spectrum is insensitive to these anomalies, 
it is necessary to perform additional investigations of 
the temperature sky map. These are d o ne in harmonic, 
wavelet, and pixel space ( Cabella et all |2004 ). Even if 
Gaussianity holds, it may still give new insights to switch 
to another representation of the statistical properties of 
the temperature maps since a phenomenon can be more 
easily detected in one representation than in another. 

The goal of this work is to provide a clear and intuitive 
analysis in pixel space regarding abundances of large-scale 
hot and cold spots identified as regions whose mean tem- 



perature contrasts exceed some (variable) threshold. We 
analyze both observed CMB maps and Gaussian simu- 
lations based on ACDM. The comparison reveals severe 
deviations. Other authors who worked with statistics of 
local extrema in t he tempera t ure field also observed sig- 
nificant ano malies ( Hou et all , 120091; lLarson and Wandeltl . 
12004120051) . 

We start by recalling some basic results that connect 
pixel-space analyses with the angular power spectrum in 
Sec.ini A comprehensive description of our method follows 
in Sec. IIIII including the preparation of adequate Gaussian 
simulations, the working principle of our spot searching al- 
gorithm, and an error estimation. Our results are presented 
in Sec. [TV] We consider both cut-sky maps (with unreliable 
pixels excluded by the KQ75 temperature analysis mask) 
and the Internal Linear Combination (ILC) full-sky map, 
and quantify deviations from Gaussian simulations. We 
sum up and conclude in Sec. |Vl 



II. PRELIMINARY CONSIDERATIONS 

The most robust comparison between predicted and ob- 
served spot abundances of CMB sky maps relies on simu- 
lated maps since analytic methods can hardly care for com- 
plications due to masking and beam properties. Creating 
a number of simulated maps and treating them in exactly 
the same way as the original map therefore is the clear- 
est method. Nonetheless, it is instructive to recall some 
well-known analytic results that connect the pixel-space 
analysis to familiar harmonic space. 

The spot abundances in a CMB sky map are dictated by 
the angular correlations of temperature fluctuations. The 
most popular theories stick to Gaussianity and statistical 
isotropy. Then, the ensemble average of the angular cor- 
relation between two directions (6, ip) and (8', if') only de- 
pends on the angle O between them. This leads to the 
definition of the angular correlation function 



C(6) 
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We can switch to harmonic space by decomposing the tern- 



2 



perature field into spherical harmonics: 



(2) 



Lm 



where the crucial assumption of statistical isotropy implies 

(fli m a\, m ,) = 8a>5 mm >Ci. (3) 

So, in this case, all the statistical information is in the co- 
efficients Cg, the angular power spectrum. More generally, 
we may define 



2f 



(4) 



When searching for spots of a given size, we will aver- 
age the temperature fluctuations in regions of that size. 
These regions are defined by window functions W(6,ip). 
The mean temperature contrast in such a region is 



AT = J dn AT(9,<p)W(6,<p). 



(5) 



In our sense, a spot is characterized as follows. When a 
threshold AT is fixed, a hot spot is found if AT > AT, 
whereas a cold spot is found if AT < —AT. The charac- 
teristic scale for AT is the mean temperature contrast for 
these regions AT rms = VlATSy. Clearly, if AT « AT rms , 
most regions will be spots, if AT 3> AT rms , only a few or 
none. 

The transformation to harmonic space can be done by 
decomposing the window function W{0, tp) into spherical 
harmonics with coefficients Wi m and defining 



(6) 



Together with Eqs. j2|) and (jSJ, it is straightforward to 
calculate 



AT 2 = V 

rms / j 



2£+l 
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(7) 



This result shows that the mean temperature fluctuation 
AT rms is given by the multipoles Ce weighted by We- The 
We strongly depend on the angular scale of the regions. 
Their magnitude will suppress large £ values correspond- 
ing to scales smaller than the window. By virtue of the 
addition theorem for spherical harmonics, we can write 

We=[dQ / dfi' W(6,<p) W(6' ,cp') P e (cosQ). (8) 



This allows us to calculate the We for a chosen window. 
An example is shown in Fig. [TJ 

In our case, it is adequate to approximate the sphere by 
the tangent plane at a region, replacing the direction (#, tp) 
by points x on the plane. For our purposes, it is most con- 
venient to work with top hat windows because they have 
clear boundaries. This is the easiest way to avoid ambi- 
guities arising from overlapping spots. Exemplary choices 
may be the top hat circle with window function 



W(x) 



1 



irR 2 



e(-R-W) 



(9) 
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FIG. 1 Coefficients We for the top hat circle window function 
at scales a = 1° (right plot), 6° (left plot). The plots show 
which multipoles predominantly determine AT rms . For smaller 
angular scale a, higher I values enter the analysis. 



or a square with window function 



1 



W(x) = -=■ 0(a - a*) Q(a - x 2 ). 



(10) 



Following iDurrerl ( 20081 . p. 218), we can approximate the 
We by an angular average over the Fourier transform of 
W(x) which considerably reduces the computational effort: 



1 f 2n 
W e ~— da\W(l)f 
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(11) 



For the aforementioned window functions, we can use this 
equation to easily calculate AT rms by Eq. ([7]). The results 
are plotted for the ACDM best-fit power spectrum in Fig. [2] 
For the sake of comparability, we use the parameter a which 
equals the square root of the windows' area; in the case of 
squares, it simply is the side length. We also show the 
relative deviation due to the different window functions. 
We conclude that the result is not sensitive to the exact 
geometry if the covered surface area is the same. 



III. METHOD 

Our strategy consists of performing an identical analysis 
of spot abundances both for observational maps and maps 
generated from simulations of Gaussian fluctuations. For 
the simulated maps, we use the best-fit ACDM model and a 
Gaussian fluctuation model based on the C'e quoted by the 
WMAP collaboration. The comparison with maps from 
observation tests Gaussianity. 

Because of the excellent data products of the WMAP 
team available at the leg acy archive 1 and the comprehen- 
sive HEALPix package 2 ( Gorski et al.1 . 120051 ). it is possi- 
ble to obtain reliable CMB sky maps and to create maps 
from Gaussian simulations. We summarize the steps in 



1 http://cmbdata.gsfc.nasa.gov 

2 http://healpix.jpl.nasa.gov 
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FIG. 2 Mean temperature fluctuation for various spot sizes and 
the ACDM power spectrum. The plots for circles and squares 
are visually indistinguishable. The difference between the result 
for circles and the result for squares is shown in the second 
figure. 



FIG. 3 The foreground-reduced V map (temperature contrast 
in mK) and the KQ75 mask cutting out the contaminated 
galaxy region and point sources. 



Sec. IIII.AI We developed an algorithm searching for hot 
and cold spots (in the sense of SecQI]) within these temper- 
ature sky maps. Its working principle and properties are 
presented in Sec. IIII.B1 The treatment of statistical errors 
is described in Sec. IIII.CI 



A. Maps and data preparation 

Whenever the original signal is to be extracted from 
CMB data, it is crucial to minimize the influence of 
foreground contamination. The frequency dependence of 
the foreground components (e.g., synchrotron emission, 
free-free emission, and thermal dust) allows to reduce 
the cont amination with t he help of various foreground 
models (|Gold et all l2009h . The WMAP team provides 
foreground-reduced maps for the Q (35-46 GHz), V (53-69 
GHz), and W (82-106 GHz) bands. Since the V band has 
a better signal-to-noise ratio than the W band and is less 
foreg round contaminated than the Q band ( Hinshaw et all 
1200711 . it is the natural choice to use the foreground-reduced 
V map. Further noise minimization by constructing lin- 
ear combinations of the maps is possible but does not af- 
fect our analysis which focuses on large scales. But still, 
large parts of the temperature map are unreliable and must 
be excluded from the analysis. We therefore apply the 
KQ75 mask, cutting out the contaminated galaxy region 



Gaussian simulations based on some input Ci spectrum 
and a beam window function are achieved with the help 
of the SYNFAST HEALPix facility. These input data can 
be obtained from the legacy archive. The power spectra 
we used are the ACDM best fit and the original W MAP-5 
estimate both shown in Fig. 1 o fNolta et al. (2009). Subse- 
quently, we will refer to them by "ACDM" and "WMAP-5" 
power spectrum for short. We take care of treating simu- 
lated and original maps as equally as possible. This neces- 
sitates the additional simulation of the instruments' noise, 
masking, and removal of monopole and dipole. Since the 
WMAP design minimiz es noise correlatio n between neigh- 
boring pixels in a map ( Page et all l2003ll , it is legitimate 
to add white noise with the properties described by the 
WMAP team at the legacy archive. 

When studying possible anisotropy of the CMB, we need 
a full-sky (unmasked) map. Since the foreground contam- 
inations usually force us to mask parts of the sky, it is not 
a trivial task to reconstruct the full-sky CMB signal. How- 
ever, the WMAP team tries to tackle this job by combining 
the measurements of all ba nds and merge the m into a single 
(ILC) map of the full sky (|Gold et all l2009jl . The applied 
procedure is independent of foreground models but has the 
disadvantage of being doubtful on scales below approxi- 
mately 10° according to the WMAP product description 
at the legacy archive. But since we are lacking any bet- 
ter alternative, we employ the 5-year WMAP ILC map for 
full-sky analyses. 
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B. Spot searching algorithm 

The primary goal of the algorithm is to count hot and 
cold spots in CMB sky maps on various scales and tem- 
perature contrasts. A typical application will be to plot 
spot abundances against the threshold on the temperature 
contrast AT for a specific angular scale. This application 
directly imposes several features the algorithm should have: 

(i) It must define sectors on the sphere of equal surface 
area (for some desired scale). Their mean tempera- 
ture contrasts will decide whether they are counted 
as spots. 

(ii) The areas must be chosen such that one can smoothly 
scan through the map. Between two distinct areas, 
there must exist many others allowing for a smooth 
transition. 

(iii) Double counting of spots has to be excluded. The 
easiest way to achieve this is working with top hat 
windows which have clear boundaries. Overlapping 
spots will be counted as a single. 

(iv) For a statistically satisfactory comparison between 
observed and simulated CMB maps, the algorithm 
will have to analyze many sky maps. Given the huge 
amount of data, one has to implement the algorithm 
carefully in order to make this numerically tractable. 

The algorithm is designed such that it allows for an ap- 
proximate pixelization of the sphere into distinct areas of 
a given scale. Calculating their temperature contrasts de- 
termines the mean temperature fluctuation AT rms on that 
scale. By virtue of the ergodic theorem, this is a good 
estimate for the ensemble average introduced in Sec. HU 

1. Working principle 

The first task is to define the sectors S of equal surface 
area on the sphere satisfying the requirements explained 
above. We choose them to be intersections of latitude and 
longitude rings. A latitude ring lZ lat consists of all points 
between two latitude angles 9q and 9±, a longitude ring 
7Z lon of all points between two longitude angles ipo and 
tpi. The rings have two nice properties. First, as needed 
for spot searching, one can smoothly go from one ring to 
any other ring by smoothly changing its boundary angles; 
second, as needed for calculating AT rms , it is an easy task 
to discretize a sphere into distinct rings. Since sectors are 
intersections S = lZ lat n 7l lon , they share these properties. 
We thereby satisfy the requirement of smooth scanning to 
all directions. 

We impose [meeting the requirement (i) above] equal 
area A for all sectors: 

A= dQ = dp d9 smd. (12) 

Once we have decided to define sectors like this, we still 
have some freedom to choose the boundaries #0j fo, 
ipi. In order to avoid the influence of small scales, we 
must reasonably choose the sectors such that they are by 



no means degenerated. We therefore fix this freedom by 
adding another constraint. For any sector <S, the boundary 
lines in the north-south direction and the longer boundary 
in the east-west direction are chosen to be of equal length: 

(<Pi ~ <Po) sin0* = (0i - 9 ). (13) 

On the northern hemisphere 9* — 9\ , on the southern hemi- 
sphere 9* — 0n. Note that these sectors behave well. In the 
limiting case near the equator, they correspond to squares 
in flat space. At the poles (0n = or 9\ = tt), they become 
equilateral triangles. 

In practice though, the temperature field is not given 
as a smooth function of 9 and tp. The WMAP tempera- 
ture sky maps are lists assigning a temperature contrast 
ATi to each HEALPix pixel pi. The mapping pi (9, (p) 
is given in the form of a table. But since our approach 
defines sectors by means of angles, we need the reverse. 
Given the list pi i— >• (9,<p), finding the appropriate pixel 
Pi for given angles (9, ip) corresponds to searching through 
the list. Whereas searching in an unsorted list is very ex- 
pensive, an adequate sorting may considerably reduce the 
effort. The algorithm performs the following steps starting 
at the north pole 9q — 0: 

1. For given 9q and area A, calculate 9\ and Aip by 
solving Eqs. $FZ\) and (fISj) . 

2. Collect the pixels {pi} belonging to the latitude ring 
lZ lat between 9q and 9\ . This can be done efficiently 
if the map was prepared by transforming to sorted 
latitude angles (HEALPix ring ordering). 

3. Using a fast routine, sort the list {p{\ with respect 
to longitude angles. This new sorting allows one to 
directly identify the pixels out of {p{\ belonging to 
a longitude ring lZ lon with boundaries ipo and ipi — 
these pixels form the sector S = K tat H K ton . Start 
at tpa = and <p\ — At/? and smoothly scan (by 
increasing po, p>\ by a small step size h) through 
all longitude rings. For every sector, calculate the 
sector's mean temperature contrast AT by averaging 
over the pixel values AX^ and compare it with the 
threshold AT. If it exceeds the threshold, count a 
spot if the sector does not overlap with a previously 
found spot. 

4. Choose the next ring by slightly increasing 9q H> 
9q + h. It is profitable to exploit the fact that the 
sorting for longitude angles (point 3) need not be re- 
peated completely. The algorithm saves the previous 
sorting and uses it for a presort such that as much 
information is transferred as possible. 

Having increased the threshold AT, again searching for 
spot abundances in a map can be optimized by noticing 
that spots at a higher threshold cannot be found where 
there were not spots at a lower threshold. Our algorithm 
can focus on areas around previously found spots once this 
becomes advantageous. 

If we slightly adapt the algorithm, we can use it to mea- 
sure AT rms . Now, the algorithm jumps between distinct 
sectors instead of smoothly transforming them. The dis- 
tinct sectors are visualized in Fig. 0] In every distinct sec- 
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of equal area now concerns the new shape and reads 



d9 sm6 Wi(9,<p) = A. 



(14) 



As an example, we compare the standard shape with top 
hat circles of equal area [cf. Eq. ©] and plot the result 
in Fig. [5] For low thresholds, the abundances are system- 



FIG. 4 Exemplary decomposition of the sky into N aec distinct 
sectors Sj for measuring AT rma . For searching spots, the algo- 
rithm analyzes many more sectors S (those in between, sharing 
pixels with the illustrated sectors Sj). Nonetheless, N aec limits 
the maximum number of spots since overlapping spots are not 
multiply counted. 



tor, the mean temperature fluctuation is calculated. The 
squares are averaged to give AT rms . Although the shapes 
of the sectors vary, the results of SecHJensure that AT rms 
is only marginally affected. 
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2. Treatment of masked maps 

The sectors defined by our algorithm may include none, 
some, or many masked pixels. We must define selection 
rules determining which sectors are to be included in the 
statistics. We used the following two rules. The most re- 
strictive choice is to only consider sectors with no mask 
overlap (strict selection for short). These sectors will only 
contain reliable pixels. But since especially on large scales, 
only a minority of sectors will belong to this group, bad 
statistics are the price to pay. The alternative choice is to 
also consider sectors with a slight mask overlap (tolerant 
selection). This is a compromise between good statistics 
on the one hand and reliable results on the other. We 
typically allow for 5% masked area within a sector which 
guarantees that usually the majority of sectors fall into this 
group. In any case, we emphasize that masked pixels, even 
if included in the statistics, are assigned zero temperature 
fluctuation. This will avoid misinterpreting foregrounds 
as a CMB signal. Note however, that the pixels of zero 
temperature fluctuation reduce AT rms . For comparisons 
between observed maps and Gaussian simulations, we em- 
ploy the tolerant selection for the sake of better statistics; 
the comparison is still trustworthy. 



3. Alternative shapes 

The algorithm works with the shapes defined in 
Sec. IIII. B. II and illustrated in Fig. 2J But we can easily 
treat other shapes by embedding them into the previous 
sectors. This corresponds to a multiplication of the pre- 
vious window function Wo with the window function W\ 
of the desired shape where Wq must be large enough to 
ensure Wo = 1 where W\ is non-zero. The condition (fT2]) 



FIG. 5 Mean spot abundances in 100 simulated ACDM full-sky 
maps showing the results for different window functions of scale 
a = VA = 6°. 

atically higher for the standard window function. This is 
due to the fact that circles do not exhaust the area with- 
out space in between. The effect becomes important where 
many spots are found and overlapping is frequent but dis- 
appears for large thresholds where the results agree. 

4. Step size dependence 

In the ideal case, the boundary angles of the sectors 
would vary in a perfectly smooth manner when searching 
for spots in a map. But numerically, we have to choose 
a finite step size h (introduced in Sec. IIII.B. 1|) . A good 
choice of h balances sensitivity and numerical effort. Fig- 
ure [6] shows detected spot abundances against h in simu- 
lated maps. We chose h = 0.3° for which we conclude that 
our sensitivity to detect spots is satisfactory. 

C. Errors and cosmic variance 

There are statistical uncertainties simply due to the fi- 
nite number of Gaussian simulations. Moreover, the CMB 
signal itself can be regarded as the outcome of a statisti- 
cal process. It is therefore subject to statistical variation, 
quantified by the concept of cosmic variance. 

Let us assume that N Gaussian maps are analyzed for 
spots (area and threshold fixed). If n( fc ) spots are detected 
in map k, the mean spot abundance is 
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FIG. 6 Mean spot abundances for a fixed threshold (80 uK) 
against a varying step size. fOO masked ACDM simulated maps 
were scanned for a = 6°, the error bars quantify the statistical 
error. 



The statistical uncertainty of the mean value n and the 
statistical deviation of the single values riS k ' are 



££i(» W 



nf 



N(N — 1) 



= Nai. 



(16) 



The same procedure applies if we measure the mean tem- 

(k) 

perature fluctuations AT rmi5 in the maps and calculate a 
mean value AT^ S . 

We now consider cosmic variance. When we observe a 
spot abundance n, we must expect a certain deviation from 
the theoretically predicted ensemble average (n). The ex- 
pectation value of this deviation, cr 2 = ((rt — (rt)) 2 ), quan- 
tifies cosmic variance. For a very large number N of simu- 
lated maps, we may replace the ensemble average by an 
averaging over the set of simulations. We then obtain 
a n « c n <fc) with the latter calculated according to Eq. (TTfJl) . 
This can be done equally for the mean temperature con- 
trast AT rms . Whenever we specify cosmic variance (e.g., 
in the form of error bars), we estimated it by this method. 



A. Cut-sky maps 

The spots' size is characterized by their area A. We 
use the parameter a = yl to specify the angular scale of 
this size. Since on the one hand, we aim at large scales, 
and on the other hand, we want reasonable statistics, we 
are forced to find a compromise. We chose an angular 
scale a = 6°. The spot abundances of the WMAP-5 V 
map and 500 Gaussian ACDM simulations (created as de- 
scribed in Sec. IIII. A|) are found for varying threshold AT. 
The HEALPix resolution parameter of the maps is 8, corre- 
sponding to N P i X — 12 x 256 2 = 786,432 pixels. Statistical 
uncertainties and cosmic variance are displayed as error 
bars even though the spot abundances for different thresh- 
olds are of course correlated. The results for hot and cold 
spots are plotted in Fig. [7] The striking feature of the plots 
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IV. RESULTS 

The application of the spot-searching algorithm de- 
scribed in Sec. IIIII shows that the standard model ACDM 
together with Gaussianity predicts more large-scale hot 
and cold spots than are actually present in cut-sky 
WMAP-5 data (see Sec. IIV.A| . Removing the quadrupole 
or using the original WMAP-5 C e (instead of the ACDM 
fit) considerably reduces the discrepancies. While only 
0.16%-0.62% of Gaussian ACDM simulations fall below the 
observed mean temperature fluctuations on angular scales 
of 4°-8°, this increases to 2.5%-8% if the quadrupole is 
removed. We also investigate full-sky maps in Sec. IIV.BI 
and modifications of the first multipoles in Sec. IIV.CI 



FIG. 7 Spot abundances in the CMB sky (with cosmic vari- 
ance) as compared to 500 ACDM simulations (with statistical 
errors) on an angular scale of a = 6°. The fractions of Gaussian 
simulations with smaller values of s [Eq. (|17[) ] are p^ ot = 0.2% 
and pf M = 1.8%. 

is the discrepancy between theory and observation. They 
only agree in the limit of very small thresholds AT where 
it is obvious that almost every area is counted as a spot 
anyway. The discrepancy is seemingly more drastic for hot 
spots. In the plot for cold spots, it is seen that there is one 
considerable cold spot nearly reaching 150 \xK. But even 
this spot does not surpass the ACDM prediction. We note 
that this spot is localized in the region of the famous Vielva 
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cold spot (|Vielva et all 12004) . It is insightful to look at the 
spot distributions of single Gaussian simulations in order to 
get an impression of their typical behavior. Five examples 
are plotted in Fig. [5J 
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FIG. 8 Spot abundances in five randomly chosen Gaussian sim- 
ulations based on the ACDM best-fit power spectrum and the 
mean curve from Fig. [7J (hot spots) . 



Because of the strong correlation between the spot abun- 
dances Hi for different thresholds AXi, it is difficult to judge 
the significance of the discrepancies by eye. A possible 
quantity that can be used for a comparison of the observed 
CMB map with Gaussian simulations is obtained by sum- 
ming up the spot abundances at different thresholds, 



(17) 



TABLE I The fraction p of Gaussian ACDM simulations with 
a AT rms smaller than found in the V map on the angular scale 
a. 
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Fraction p 


4° 


0.50% 
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0.36% 



ACDM simulations (again with tolerant selection). We see 
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FIG. 9 The mean temperature fluctuation for different angular 
scales a in 50 Gaussian ACDM simulations and the VI map. 



where the lowest threshold included is chosen to be the 
characteristic scale Af rms . We denote the fraction of 
Gaussian simulations k with smaller than found in the 

V map by p s . For the spot abundances shown in Fig. [71 
we find p^ ot = 0.2% for hot spots and p c s old = 1.8% for cold 
spots. 

The discrepancies are reflected in the mean temperature 
fluctuation AT rms which on large scales is higher in ACDM 
simulations than in the observed CMB sky. We have simu- 
lated 5000 ACDM maps and compared their mean temper- 
ature fluctuations to the value of the V map. We employed 
the tolerant selection of our algorithm (see Sec. IIII.B.2I) . 
For a — 6°, we find the mean value AT rms = 39.4 pK 
for the V map, as compared to the mean value AT rms — 
47.9 ±0.1 pK for ACDM, where the error is only statistical 
while cosmic variance amounts to 4.2 pK. Only a fraction 
p = 0.6% of the simulations had a smaller AT rms than the 

V map. This does not improve at other large angular scales 
which can be seen in Table Q] It is interesting to see how 
this behavior changes when going to smaller scales. How- 
ever, the results on smaller scales (approaching 1°) become 
sensitive to noise and beam properties. Since the WMAP 
team off ers the latter fo r the differencing assemblies VI 
and V2 dHill et all . 120091) instead of the combined V map, 
it is the easiest to switch to the Fl-band map and simula- 
tions thereof. The impact on large scales is negligible. Fig- 
ure [HI shows AT rms against the scale a for the VI map and 



that the deviations decrease when going to smaller scales. 
This is also suggested by the Ct spectrum which is in good 
agreement with the ACD M fit for large I which do minate 
on small scales. But still, iMonteserin et al I (|2008l) find a 
too low CMB variance which approximately corresponds 
to AT rms on scales even smaller than 1°. 

For the results in Fig. [5J we used the highest available 
HEALPix resolution 9 corresponding to N p i x = 12 x 512 2 = 
3, 145, 728 pixels in a map. As stated above, the plots are 
highly influenced by the beam function and noise. The 
beam function acts as an extra window function which 
suppresses the growth of AT rms for small scales. The 
white noise instead leads to a diverging 1/a behavior on 
the smallest scales (with an effective pixel noise ampli- 
tude Opix and the number of pixels N a = N p i X x a 2 /47r 
within a sector of scale a, the noise contribution will be 
AT?± e = cWVAa oc 1/a). 



On large scales, the first multipoles of the Ct spectrum 
play an important role (see, e.g., Fig. [TJ. It is there- 
fore a nat ural idea to su s pect t he well-known quadrupole 
anomaly ( Hinsha w et all [2007) to be responsible for the 



observed discrepancies. We check this by repeating the 
analysis after removing the quadrupole from the ACDM 
simulations as well as the observed CMB map. The re- 
sults, summarized in Table ITl) confirm the influence of the 
quadrupole anomaly. Now, the fractions p of Gaussian 
ACDM simulations reach p = 7.3% for a = 6°. These 
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TABLE II The fraction p of 1000 Gaussian ACDM simulations 
with a AT rms smaller than found in the V map on the angular 
scale a, after removing the quadrupole from the maps. 



Scale a 


Fraction p 


4° 


6.5% 


5° 


8.0% 


6° 


7.3% 


7° 


2.5% 


8° 


6.4% 



TABLE III The fraction p of 1000 Gaussian simulations 
(WMAP-5 Ct) with a AT rms smaller than found in the V map 
on the angular scale a. 



Scale a 


Fraction p 


4° 


4.2% 


5° 


5.4% 


6° 


5.8% 


7° 


2.4% 


8° 


4.1% 



numbers still do not show good agreement, but they are 
not statistically significant anymore. 

We now investigate whether there are still discrepancies 
if we compare the observed V map with Gaussian simula- 
tions based on the original WMAP-5 Cg spectrum rather 
than the ACDM best fit. This tests whether the observed 
map is a typical Gaussian realization of the WMAP-5 
power spectrum. Figure [10] shows the spot abundances. 
The effect arising from changing the power spectrum is 
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FIG. 10 Spot abundances in the CMB sky (with cosmic vari- 
ance) as compared to 500 simulations (with statistical errors) 
based on the WMAP-5 Ce spectrum on an angular scale of 
a = 6°. The fractions of Gaussian simulations with smaller 
values of s [Eq. (JT7J)] are p* ot = 3.4% and p c ° M = 13.2%. 



clearly visible and reduces the discrepancies to some ex- 
tent. But although closer to the spot abundances in the 



observed cut-sky CMB map, the numbers of hot and cold 
spots arc still too high. Again, this is reflected in the fact 
that most simulated maps have a larger AT rms than the V 
map. Even though the values, listed in Table IHI[ are less 
drastic, we emphasize that the WMAP-5 estimation of the 
Ce relies on similar data, i.e., cut-sky CMB maps. If the 
observed CMB map was a typical Gaussian realization of 
the extracted Cg spectrum, we would expect agreement. 

Bearing in mind, however, that power spectra refer to 
the full sky whereas we only look at regions outside the 
mask, an explanation could be that the missing spots were 
located in the hidden part of the sky. In the next section, 
we investigate whether the WMAP-5 ILC map indicates 
this violation of isotropy. 



B. ILC full-sky map 

The five-year ILC map is the best approximate full-sky 
CMB map available. We therefore analyze it even though 
the quality of the reconstruction is not high enough to guar- 
antee robustness of the results (see, also, Sec. IIII.A"|) . We 
analyze the ILC full-sky map and 100 Gaussian full-sky 
simulations based on the WMAP-5 power spectrum and 
separately consider the results in three sky regions. First, 
we analyze the full sky. Second, we collect the spots of 
those regions that have also been studied in the V map, i.e., 
regions with no or little overlap with the KQ75 mask (tol- 
erant selection). Finally, we consider the remaining spots 
that consequently lie in sectors completely inside the mask 
or with considerable mask overlap (rejected by tolerant se- 
lection). We loosely refer to the three regions as full sky, 
outside, and inside mask. The results are plotted in Fig.lTTl 
In the previously analyzed region (outside the mask), we 
see too few spots, as before. But there are by far too many 
spots in the complementary region. The variances provid- 
ing the error bars are, as always, obtained from Eq. (|16[), 
Although there are less statistics inside the mask than in 
the full sky, the error bars in the corresponding figure are 
smaller. This can be intuitively understood as follows. If, 
for simplicity, we assumed that the spot abundances out- 
side and inside the mask were statistically independent, the 
variances of n , a^ ut would add to of ull in the full sky, whence 
fin < <7fuU- The loss of statistics when counting spots in- 
side the masked region only causes the relative fluctuations 
between two Gaussian simulations to increase. The error 
bars in the central figure (outside the mask) visually ap- 
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conclusion (iBernui and Reboucasl I2009I : ICopi et all . 2009: 
iHaiianl 1200/1) . 

Anisotropy of the CMB is a possible explanation of the 
discrepancies revealed in Sec. IIV.AI and quantified in Ta- 
ble HI and indeed, the ILC map contains this anisotropy. 
But since there is not enough reliable information about 
the CMB signal in the galactic plane, we cannot finally 
judge whether this is the true solution to the problem. We 
have also studied if, additional to the galactic plane, the 
orientation of the galactic halo defines a preferred direc- 
tion. Therefore, we divided the ILC map into two halves, 
one around the galactic center and one covering the oppo- 
site direction. We have seen no signal of anisotropy in this 
direction. 
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FIG. 11 Spot abundances in the ILC map (with cosmic vari- 
ance) compared to simulations (with statistical errors) based on 
the measured WMAP-5 Ci on an angular scale of a = 6° for 
three different parts of the sky. The corresponding values of p s 
are p{ uU = 58%, p° ut = 7%, and pf = 96%. 



pear larger due to the logarithmic plotting but are in fact 
smaller than for the full sky. 

The values of p s confirm the uneven distribution of 
spots in the ILC map. For the full sky, we have p* u11 = 
58% in good agreement with the simulations. Outside 
the mask, there are too few spots, p° ut — 7%, whereas 
inside the mask, we find p' s n = 96%. The ILC map 
is clearly anisotropic. Other authors draw the same 



C. Modified power spectra 

We have pointed out that anisotropy is a potential expla- 
nation. It is however unsatisfying to assume that so many 
additional spots lie in the contaminated regions hidden by 
the KQ75 mask. This would be a surprising coincidence of 
CMB signals and the orientation of the galactic plane. Al- 
ternatively, we may stick to statistical isotropy; then, our 
results may be due to some non-Gaussian signal. 

In this section, we investigate whether our results imply 
non-Gaussianity or statistical anisotropy by themselves. 
We do this by analyzing the effect of modifications to the 
Ci spectrum. 

So far, AT rms has proved to be a good parameter to 
quantify the visible effects. We can perform a quick check 
whether our data supports the hypothesis that AT rms is 
the decisive parameter. Out of the 500 simulations with 
the WMAP-5 power spectrum used in Sec. IIV.AI we collect 
those with a AT rms smaller or equal than those found in 
the V map. Figure [T2l shows their spot abundances which 
agree well with the V map. 

If there is a Cg spectrum that produces AT rms values 
similar to the ones found in the V map, our results alone 
do not imply non-Gaussianity or statistical anisotropy. In 
order to keep the analysis as generic as possible, we do 
not use any specific cosmological model but only modify 
the Ci of the original WMAP-5 spectrum. Figure [§] sug- 
gests that only large scales are affecte d which is why we 
concentrate on a few low multipoles I. ICopi et all (|2009T) 
found that the correlation function is essentially zero on 
angular scales above ~ 60°. Since this scale is roughly 
linked to the multipole range t < 3, our first modifica- 
tion simply consists in setting Ci = for £ < 3 (although 
of course the correlation function does not translate this 
easily). Another example may be to halve the Ci for 
I < 5 (modification II) . Figure Q~3] shows the resulting val- 
ues of AT rms . The plots show the discrepancies between 
the ACDM prediction, the WMAP-5 spectrum, and obser- 
vation. We also show the results for a combined power 
spectrum, replacing the first 32 multipole s by the values 
quoted by WMAP-1 (jHinshaw et al.l . l2003l) . For this range 
of multipoles the WMAP analysis changed aft e r the 1-year 
release, following the suggestion of lEfstathioul ( 2004) . The 
difference between WMAP-5 and WMAP-1 may serve as 
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FIG. 12 Spot abundances of Gaussian simulations k (errors 
statistical) with AT^ S < ATX^* P in comparison to the V map 
(with cosmic variance). For these plots, we have p^ ot — 68% 
and Ps° ld = 86%, showing agreement. Since we only consider 
Gaussian simulations with ATrms smaller than in the V map, 
it is no surprise that the p s values lie above 50%. 



an illustration that the extraction of reliable Ce for low i is 
a nontrivial task. Modifications I and II of the power spec- 
trum succeeded in reconciling Gaussian simulations and 
observed CMB sky. This is confirmed by measuring the 
spot abundances in simulated maps based on the modified 
spectra, seen in Fig. Q3] We conclude that our results are 
not incompatible with Gaussianity. However, if we stick to 
Gaussianity, they favor (although statistically not signifi- 
cant, cf. Table ITU)) even lower values of the first multipoles 
than currently estimated. 



V. CONCLUSIONS 

The study of spot abundances has revealed discrepan- 
cies between the cut-sky CMB temperature maps and the 
standard best-fit ACDM model or, but less significant, a 
Gaussian spectrum for the Ce estimated by WMAP-5. We 
have shown in Sec. lIV.Cl that a good parameter to quantify 
them is the mean temperature fluctuation AT rms which we 
investigated on large scales. On scales a between 4° and 
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FIG. 13 The mean temperature fluctuation for large angular 
scales a. We compare the VI map with ACDM simulations 
(highest AT rrns ), simulations based on the WMAP-5 Ce (a bit 
lower), the WMAP-1 Ce for I < 33 (still lower), and on two 
modified spectra. The first modification is created by setting 
Ce — for I < 3, the second by halving the Ce for £ < 5. The 
modified spectra agree well with the VI map. 



8°, only 0.16% to 0.62% of Gaussian simulations based on 
the ACDM best-fit power spectrum fall below the AT rms 
value of the observed V map. If this merely was an imprint 
of the anomalously low quadrupole, we would expect the 
discrepancies to disappear when removing the quadrupole 
from the Gaussian simulations and the V map. The differ- 
ence in fact reduces, the aforementioned fractions change 
to 2.5% to 8.0%. These numbers are not significant and 
do not allow for a clear interpretation whether our results 
go beyond the quadrupole anomaly. Similar fractions are 
obtained when exchanging the ACDM best-fit spectrum 
by the originally published WMAP-5 Ce, yielding 2.4% to 
5.8%. This is difficult to understand if we bear in mind that 
the C f themselves are e stimated from the cut-sky CMB 
maps (|Nolta et all . 120091) . 



Non-Gaussianity and also statistical anisotropy are pos- 
sible explanations. In our case, anisotropy means that 
many spots have to be hidden behind the masked re- 
gion. Unfortunately, this hypothesis can hardly be tested 
as there is currently no method to reliably extract the 
CMB signal in the highly foreground-contaminated regions. 
Nonetheless, we have employed the WMAP-5 ILC full-sky 
CMB map and found evidence for anis otropy i n this map. 
This agrees with results obtained by iHaiianl ( 20071 ) and 
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FIG. 14 Spot abundances in the CMB sky (with cosmic variance) as compared to 100 simulations based on the two modified 
spectra, respectively, (errors statistical) on an ang ular scale of a = 6°. The first modification yields p* ot = 30% and p c ° M = 55%, 
the second modification p^ ot = 37% and p c ° ld = 67%. 



Copi et al.l (2009) who found that most of the power on 
the largest scales comes from the (masked) galaxy region. 
Though possible, this unnatural alignment of the CMB sig- 
nal with the galactic plane would be intriguing and lacks 
so far any explanation. 

Our analysis of Sec. IIV.CI shows that our results for 
cut-sky maps do not suggest non-Gaussianity or statistical 
anisotropy by themselves. They agree well with Gaussian 
fluctuations if one performs a modification of the lowest 
multipoles. In doing so, no fine-tuning of the Ce is neces- 
sary in order to reconcile the spot abundances from Gaus- 
sian simulations and the observed CMB. It is sufficient to 
lower the first multipoles by a substantial amo unt. When 
study ing local extrema in the temperature field, iHou et ahl 
(|2009h similarly found discrepancies that disappeared when 
excluding the first multipoles. We recall, however, that the 
Ce and the assumption of Gaussianity completely fix the 
expected spot abundances. If both the extraction of the Ci 
by WMAP-5 and our analysis of spot abundances are cor- 
rect, our results may indicate non-Gaussianity or statistical 
anisotropy. 

If the discrepancies are not caused by mere statisti- 
cal coincidence or unknown secondary effects, we have to 
leave open the question whether we see the consequence 
of non-Gaussianity or anisotropy, or whether our results 



strengthen the evidence for a severe lack of large-scale 
power. The first case would challenge fundamental assump- 
tions, the second would make it difficult to understand the 
CMB maps on large scales within standard ACDM cosmol- 
ogy. If the discrepancies between the CV, as determined by 
WMAP-5, and the spot abundances persist, this can be 
interpreted as a signal for non-Gaussian fluctuations. 
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